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^ _ ^ Self-interaction corrected local spin density approximation calculations were performed for 

' La(i_^-)Sr^Mn03 (LSMO) (0.0 < x < 0.5). The influence and inter-relationship of Sr doping, 

I magnetic structure, O displacements and phase segregation on the Mn charge state were studied. A 

, half-metallic state was obtained for LSMO with manganese configuration Mn^+, whilst Mn*"'" gave 

■ rise to a metallic state with a negligible spin polarisation at the Fermi level. Elongating the MnOe 

2f^' octahedron led to a static mixed valence Mn"^"'"/Mn'*"'" configuration. In the mixed valence state the 

1^ , total energy was minimized by an ordered array of Mn*"*" and Mn'^"'" Mn02 planes which showed 

charge ordered stripes. 

PACS numbers: 75.47.-m 
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I. INTRODUCTION 



The half-metallic properties of La(i_a;)Sra:Mn03 (x=0.3) (LSMO) are of great importance for applications in spin- 
^ ironies. The tunnel magnetoresistance junction of LSMO/SrTiOa/LSMO shows magnetoresistance ratio in excess 

of 1800%i. This result, according to Ref. strongly underlines the half-metallic nature of mixed-valence mangan- 
^ ites. The electronic properties of LSMO, as described by band theory, are nearly half-metallio2i2i^, reflecting the 
I \ so-called transport half-metallic behaviour. ^-^ However the fascinating electronic and magnetic properties of LSMO, 
' O ' including colossal magnetoresistance (OMR), indicate that the electronic structure is more complex than the stan- 
^ ■ dard band theory picture (see reviewsSiS). In particular, the electronic structure is determined by the competition 
of double exchange and superexchange interactions, charge/orbital ordering instabilities, and strong coupling to the 
lattice deformations. Local Jahn- Teller effects, such as random Jahn- Teller distortions of the MnOe-^° octahedra, as 
well as dynamical effectsii, have recently been invoked to explain the magnetoconductivity and optical conductivity 
respectively. Despite the numerous studies of the phase diagram of LaMnOs-SrMnOa (LMO-SMO) (such as in Refs. 
Ci 112 13.14), there are still many conflicting interpretations of the role of the Jahn- Teller effect in this materiaUa*i*, 
localisation of d electrons and polarisation of electrons at the Fermi level. 
I In this paper we discuss issues concerning the charge ordering, and more specifically, the distribution of Mn'^"^ and 
00 . Mn*"*" in LSMO. We use the first principles self-interaction corrected local spin density (SIC-LSD) approximationi2£ 
This method can determine the number of valence band states. Hence, it can differentiate between Mn contributing 
, three states (Mn'^"'") to the valence band with the remaining four Mn states (3t2g and le^) localised well below the 
' valence band, and Mn contributing four states (Mn^+) to the valence band, with the remaining three Mn states (3t2g) 
+^ ' localised below the valence band. This method was successfully applied to the study of orbital order in LaMnOsi^ 

a Calculations with the SIC-LSD as a function of Sr doping are presented in this paper for Lai_a;Sra;Mn03 where 
0.0<x<0.5. The Sr doping is modeled with supercells and also with the rigid band approach. In particular we will be 
concerned with determining the Mn valency as a function of concentration x. It is found, as expected, that the valency 
' changes as a function of Sr doping from Mn'^"'" to Mn'*"'". However, a mixed phase of Mn^+ and Mn^+ valencies is found 
O ■ to accompany the valency change. The Mn'^"'"/Mn^+ ordering in this mixed phase is consistent with charge ordered 
O ■ stripes. Furthermore it is found that, in the Mn'*+ state and as a function of Sr doping, a change of magnetic structure 
from ferromagnetic to anti-ferromagnetic takes place. To highlight the effects of Sr doping on the Mn valencies the 
lattice parameters are kept constant for all concentrations of Sr doping. However, a slightly larger lattice parameter 
makes the Mn'^"'" valency less unfavourable. The structural properties of LSMO show an equally rich variety of charge 
ordered states. Elongating the MnOg octahedron makes it for the Mn ion more likely to take on the Mn^+ valency. 

The paper is organized as follows. In the next section we introduce the theoretical background of our electronic 
structure calculations, and in particular, the self-interaction corrected local spin density method. The technical and 
computational details regarding the application of the SIC-LSD to LSMO are discussed in Section HI. In Section 
IV, we verify that the method works separately for LaMnOa and SrMnOs. Section V discusses the correlation 
between magnetic structure and charge order in LSMO. The influence of oxygen displacements on the charge order is 
investigated in Section VI. Section VII presents the results of calculations for various realizations of phase separation 
and their influence on the charge order. The conclusions of the paper are summarized in Section VIII. 
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II. THEORY 

The basis of the SIC-LSD formaUsm is a self-interaction free total energy functional, E^^"^ , obtained by subtracting 
from the LSD total energy functional, E^^^ , a spurious self-interaction of each occupied electron state tpcT^, namely 

occ. 

a 

Here a numbers the occupied states and the self-interaction correction for the state a is 

S^J^ = U[n^]+E^,^^[n^], (2) 

with [/[ria] being the Hartree energy and E^^^[na] the LSD exchange-correlation energy for the corresponding charge 
density Ua and spin density . The SIC-LSD approach can be viewed as an extension of LSD in the sense that the 
self-interaction correction is only finite for spatially localised states, while for Bloch-like single-particle states E^^'~' 
is equal to E^^^ . Thus, the LSD minimum is also a local minimum of E^^'^ . A question now arises, whether there 
exist other competitive minima, corresponding to a finite number of localised states, which could benefit from the 
self-interaction term without loosing too much of the energy associated with band formation. This is often the case 
for rather well localised electrons like the id electrons in transition metal oxides or the 4/ electrons in rare earth 
compounds. It follows from minimisation of Eq. |^ that within the SIC-LSD approach such localised electrons 
move in a different potential than the delocalized valence electrons which respond to the effective LSD potential. For 
example, in the case of manganese, three (Mn'*+) or four (Mn'^+) Mn d electrons move in the SIC potential, while all 
other electrons feel only the effective LSD potential. Thus, by including an explicit energy contribution for an electron 
to localise, the ab-initio SIC-LSD describes both localised and delocalized electrons on an equal footing, leading to a 
greatly improved description of static Coulomb correlation effects over the LSD approximation. 

In order to make the connection between valence and localisation more explicit it is useful to define the nominal 
valence as 

Nval = Z - Ncore - NglC, 

where Z is the atomic number (25 for Mn), Ncore is the number of core (and semi-core) electrons (18 for Mn), 
and Nsic is the number of localised, i.e., self-interaction corrected, states (either three or four for Mn^+ and Mn'^+ 
respectively). Thus, in this formulation the valence is equal to the integer number of electrons available for band 
formation. To find the valence we assume various atomic configurations, consisting of different numbers of localised 
states, and minimise the SIC-LSD energy functional of Eq. ^ with respect to the number of localised electrons. The 
SIC-LSD formalism is governed by the energetics due to the fact that for each orbital the SIC differentiates between 
the energy gain due to hybridisation of an orbital with the valence bands and the energy gain upon its localisation. 
Whichever wins determines if the orbital is part of the valence band or not and in this manner also leads to the 
evaluation of the valence of elements involved. The SIC depends on the choice of orbitals and its value can differ 
substantially as a result of this. Therefore, one has to be guided by the energetics in defining the most optimally 
localised orbitals to determine the absolute energy minimum of the SIC-LSD energy functional. The advantage of 
the SIC-LSD formalism is that for such systems as transition metal oxides or rare earth compounds the lowest energy 
solution will describe the situation where some single-electron states may not be Bloch-like. For Mn, these would be 
the Mn 3d states, but not the O 2p states, as trying to localise the latter is energetically unfavourable. 

In the present work the SIC-LSD approach, has been implemented^'' within the linear muffin-tin-orbital (LMTO) 
atomic sphere approximation (ASA) band structure method,^* in the tight- binding representation^^. 

III. CALCULATIONAL DETAILS 

We performed SIC-LSD calculations for both La(i_2.)Srj;Mn03 for 0.0 < a; < 0.5 and SMO. The space group 
symmetry of the LSMO structure changes, as a function of Sr concentration, from Pnma (0% of the Sr) to i? — 3c 
(between 10% and ~50% of the Sr) and PQ3/mnic for SMO. However, for the sake of comparison, we performed 
calculations using the same cubic crystal structure for all the concentrations. Thus our results for LMO refer to a 
hypothetical cubic phase. In order to highlight the effect of the electron doping by substituting Sr for La we also kept 
the lattice parameter constant at 7.32 atomic units (a.u.) - which is the average of the theoretical lattice parameters 
of ferromagnetic LMO and SMO. The experimental lattice parameter, in going from LMO to SMO, changes by 1%, 
i.e., from 7.45 to 7.38 a.u.. 
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TABLE I: Energy for different magnetic configurations for cubic SrMnOs referred to the ground state energy. The energy 
differences are per chemical unit cell and the lattice parameter was taken to be 7.32 a.u.. The eg(3^2_r2') orbital is written as 
leg. Also displayed are the magnetic moments (MM) for each of the magnetic structures. 

AF-A AF-G FM 

AE MM AE MM AE MM 

[mRy] [(ib] [mRy] [^b] [mRy] [/ib] 

LSDA 115 2.60 108 2.49 120 2.59 

SIC(3t29 ) 5 2.74 2.79 1 2 2.68 
SIC(3t29+le9) 114 3.44 116 3.47 105 3.44 

For the linear mufRn-tin basis functions, we used 6s, bp, bd partial waves for the La and 5s, 4p, 4c? for the Sr 
atoms, and treated them as low-wavesi^ Including also 4/-basis functions on the lanthanum, treated as intermediate 
waves, was of no substantial importance for the final results. On the manganese atoms only 4s and 3d partial waves 
were treated as low waves, whilst the 4p-waves were treated as intermediate. On the oxygen only 2s and 2p partial 
waves were treated as low-waves, and 3d-waves as intermediate. The atomic sphere radii were 4.0, 2.3 and 1.8 a.u. 
for the lanthanum and strontium, manganese and oxygen, respectively. These spheres were chosen to minimise the 
discontinuity in the Hartrec potentials, giving an overlap volume of approximately 8%. No empty spheres were 
used for the cubic system. Care was taken to ensure that the results were converged with respect to both the size 
of the screening cluster and the number of k points for which the one electron equations were solved. This was 
imperative to allow comparisons to be made between different magnetic structures, which entailed the use of different 
unit cells. The screening clusters consisting of 111 atoms for the lanthanum, strontium and manganese sites and 
99 atoms for the oxygen atoms were used. The number of k points used was 256 in the full Brillouin zone for the 
paramagnetic/ferromagnetic, G-type and A- type magnetic structures. 

To perform calculations for LSMO with < x < 0.5, we have utilized the SIC-LSD within the rigid band 
model and supercell approach. In the rigid band model, the variations in the band filling (reduction by up to 0.5 
electrons) and lattice constant were the only variables depicting the change from the cubic LaMnOa to LSMO. 
Thus taking into consideration the dependence on the lattice constant, we could use the results from rigid band 
model to describe properties of other perovskites where Sr may be replaced by Ca or Ba. Supercells of the form 
La„SrmMn(„+,„)03(„+,„) (where n = 1,..,7 and m=l,2,3) were constructed to describe charge ordering effects due 
to Mn'*+ occurring in the vicinity of Sr, and Mn'^^ present around La sites. 



IV. THE END MEMBERS OF PHASE DIAGRAM: LMO AND SMO 

We start with the application of the SIC-LSD to the end compounds of LSMO, namely LaMnOa and SrMnOs. 
TableUsummarizes the LSD and SIC-LSD results for SMO for three different magnetic structures, and two Mn valence 
configurations. As seen in the table, the Mn'*+ configuration in the G-type anti-ferromagnetic (AF-G) structure is 
the ground state solution. The latter is insulating, with a band gap of 1.08 eV, twice smaller than the measured band 
gap^^ of 2.3 eV. The ground state configuration corresponds to the localization of the three t2g electrons. Localizing 
an extra d electron, the eg one with the symmetry d322_r2, giving rise to Mn'^"'', is unfavourable by more than 100 
mRy, for all three different magnetic structures. This energy difference between Mn^+ and Mn'^+ configurations in the 
ferromagnetic (FM) SMO system decreases from 120 mRy to 93 mRy, when increasing the lattice parameter from 7.2 
a.u. (corresponding to the theoretical pseudocubic FM ground state of SMO with Mn'^^) to 7.32 a.u.. On the other 
hand. Table ^ also shows that only 5 mRy separate the G-type from A- type antiferromagnetic state. This shows that 
the charge ordering energy in SMO is a much larger energy scale than the magnetic order and that its dependence 
on the lattice constant is small. The charge transfer of l.Oe and O.le respectively from Mn and O atoms to Sr atom 
hardly changes for the three different magnetic structures. For each magnetic structure the localisation increases the 
Mn magnetic moment, which, however does not change much between the different magnetic structures. 

The FM cubic LaMnOa with the lattice parameter of 7.43 a.u. has the Mn'^"'' configuration: 15 mRy separate this 
ground state from the Mn^+ excited state. In comparison, for the Jahn- Teller distorted LaMnOa structure, we find 
20 mRy energy difference i>2l Reducing the lattice constant by 1.5% to 7.32 a.u., we find that LaMnOa becomes nearly 
tetravalent and less than 5 mRy separate the Mn'^+ ground state from the Mn'*+ excited state. We find the crossover 
between trivalent and tetravalent manganese at a volume of elementary cell equal to 215 A'^. This is close to the 
expected volume of 210 from pressure experiments-^ . The charge transfer from Mn and O atoms (respectively 
0.9e and 0.2e) gives 1.5 more electrons on the La atom. Thus the charge transfer is similar between SMO and LMO. 
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TABLE II: Energy for three different magnetic configurations with reference to ground state energy for cubic Lao.sSro.sMnOs. 
Energy differences are per chemical unit cell and the lattice parameter was taken to be 7.32 a.u.. AEjibm refers to the energy 
difference from a rigid band model (RBM), whilst AEsc is the energy difference obtained using a supercell (SC) model. The 
symbol MM stands for magnetic moment. Four different valence configurations are considered: LSDA, Mn*"*" with the three 
t2g electrons localized and the two Mn'^'*' valence configurations, localizing each of the eg separately. 





AEhbm 


MM AEsc MM 




[mRy] 




[mRy] 


[/is] 


AF-A 


LSDA 


113 


2.67 


108 


2.77 


SIC(3t29) 





3.04 





2.99 


SIC(3t2s+eg(3^2 


_^2)) 63 


3.51 


63 


3.51 


SIC(3t29+eg(^2_ 


.,2)) 48 


3.53 


47 


3.53 


AF-G 


LSDA 


120 


2.62 


122 


2.56 


SIC(3t2g) 


18 


3.02 


22 


2.89 


SIC(3i29+eg(3^2 


-r2)) 64 


3.50 


62 


3.50 


SIC(3t29+e3(^2_ 


.,2)) 64 


3.50 


63 


3.50 


FM 


LSDA 


111 


2.69 


113 


2.70 


SlC{3t2g) 


4 


2.91 


6 


2.86 


SIC(3t2g+eg(3^2 


_^2)) 56 


3.51 


57 


3.50 


SIC(3t29+eg(^2_ 


.,2)) 52 


3.52 


54 


3.52 



The calculations show that the energy scales between the Mn'^+ and Mn^+ charge ordered state are much smaller 
for LMO than for SMO, specifically 20 times smaller. Actually, this is an energy scale comparable to the energy 
difference between magnetic structures. The competition between these energy scales in the case of Sr doping will be 
the subject of the next three Sections. 

V. MAGNETIC STRUCTURE AND CHARGE ORDER IN LSMO 

In Fig^we present the phase diagram of anti-ferromagnetically and ferromagneticUy ordered LSMO for Sr concen- 
tration between 0.0 and 0.5. For LMO, we obtain a ground state of Mn'^^ valency in an AF-A magnetic structure. In 
the range of x up to approximately 0.2 (excluding 0.0), the FM supercell with mixed valence of Mn'^+/Mn''^+/Mn'^+ 
gives the state with the lowest energy (see also inset of Fig^J. A crossover, as a function of Sr doping, from a FM 
mixed valence Mn'^+/Mn^+/Mn'^+ ground state to a FM Mn'^^ ground state occurs around 20% Sr doping. For Sr 
concentrations between 20% and 35% the ground state has the valency of Mn^+ in a ferromagnetic structure. For 
Sr concentrations larger than 40%, the magnetic structure changes to AF-A, but the valency remains Mn^+. For x 
larger than 0.2, we find that the total energy scales linearly between Mn'*+ and Mn'^+ configurations. We found this 
by using supercells of up to 8 formula units which allowed us to model different distributions of Mn"'+ and Mn''+ 
atoms. 

From the top part of Fig^ we note that for the unfavourable Mn'^"'' configuration the energy separation between the 
two anti-ferromagnetic structures is small. In the Mn"*^ configuration, the energy separating these anti-ferromagnetic 
structures becomes larger. Obviously, three Mn band electrons are insensitive to the magnetic structure whilst four 
band electrons, including all Cg electrons, probe the magnetic structure. From the bottom part of FigQ] we note that 
the Mn'*+/Mn'^+/Mn^+ and the Mn^+ZMn'^^ configurations are close in energy and also close to the ground state. 
Increasing the Mn'^+ concentration as realized in the Mn'^+/Mn'*+/Mn'^+ configuration takes the total energy further 
away from the ground state energy as the hole doping increases through the increase of Sr concentration. Also for 
a slightly larger lattice constant (7.43 a.u.) the deviation from the ground state energy steadily increases upon hole 
doping. The Mn"^+ configuration is energetically the most unfavourable. The differences between supercell description 
and rigid band are small on the energy scale of the separation to the Mn''+ ground state and the larger the hole doping 
becomes the more unfavourable the Mn'^"'" configuration turns out to be. 

The Lao.sSro.sMnOs system was studied using both the rigid band model and LaSrMn206 supercell. From Table 
im we note that in the LSD we obtain a FM state in the rigid band model, whilst an AF-A state is obtained in the 
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supercell approximation. The energy differences between these magnetic structures are 2 and 5 mRy for the rigid 
band model and superceU, respectively. The SIC-LSD calculations give the same ground state both in the rigid band 
model and supercell. Similarly to LSD, small energy differences of 4 and 6 mRy separate the AF-A ground state 
from the FM state in the rigid band model and supercell respectively. Figure also confirms that the rigid band 
model is an adequate description of the disordered La/Sr system. A comparison is made with supercell calculations 
of La„SrMn(„_|_i)03(„+i) with n = 1,..,7. From energy differences between Mn'*+ and Mn'^+ systems for each n 
we deduce an Mn'^+ ground state configuration for n larger than 9 (for such system as LaioSrMnnOaa) which is in 
agreement with x < 0.1 from the rigid band model. This shows the close agreement between these calculations. 

From Table ITTl we also note that the energy differences between Mn'*+ and Mn'^+ are more or less constant with 
respect to the magnetic structure and also with respect to the symmetry of the localized Cg electron in the Mn'^+ 
configuration. This energy is approximately 60 mRy. Again the energy differences between magnetic structures are 
an order of magnitude lower. 

The energies in Fig. ^ are closely balanced and the value of the lattice constant will therefore be important. 
Increasing the lattice constant by 1% to 7.39 a.u., we find that for LMO the ground state remains AF-A with Mn'^^ 
for Sr concentration x up to 0.1, and turns into FM with Mn'*+ configuration for x larger than 0.1. The Mn'^"'" 
configuration has at the same time become slightly less unfavourable. This we note, for example, by comparing the 
energy difference for Lao.ySro.sMnOa between FM Mn'*"'" and Mn'^+ configurations which reduces from 32 mRy at 
lattice constant 7.32 a.u. to 21 mRy at the 1% increased lattice constant of 7.39 a.u.. For Lao.sSro.sMnOa even the 
ground state changes from AF-A at the lattice constant of 7.32 a.u. to FM at 7.39 a.u., the energy differences being 
4 mRy and -2 mRy respectively. This shows that the combined effect of band filling and the value of the lattice 
constant gives rise to a rich variation in magnetic structure and charge order. 

The spin magnetic moments increase only slightly in comparison with LSD results. For Lao.jSro.sMnOs the Mn 
spin magnetic moment is 3.14 /x_b in SIC-LSD and 3.03 hb in LSD and for Lao.sSro.sMnOa these are 2.91 and 2.69 
/i B respectively. The density of states for LSD calculations and the SIC-LSD in Mn^'^ and Mn^+ configurations are 
shown in Fig. |21 We clearly see the majority Mn t2g peak, which occurs in the LSD calculation just below the Fermi 
level, which moves down in energy below the bottom of the valence band for the calculations for Mn^+ and Mn'^^ 
valencies. In the Mn^'^ ground state configuration for Lao.TSro.sMnOs a metallic state is obtained and the electronic 
structure in the vicinity of the Fermi level is similar to the LSD, i.e., it is also a nearly half-metallic system. However, 
for the Mn'^+ configuration, we obtain a half-metal with 1.6 eV band gap in the majority spin channel. We will see 
in the following that the occurence of a half-metallic density of states is closely associated with the Mn^+ valency. 

The importance of local effects was studied with supercells. From calculations for LaSrMn206 we note that the 
O atom in SrO layer (Osr) lost twice as many electrons as the O atom in the LaO layer (Olo)- This leads to an 
increased spin moment of 0.2 fiB for Osr in comparison with 0.07 fiB for Olq. From calculations for La2SrMn309, we 
find that Mn sandwiched between two LaO layers loses 0.07 electrons more than Mn between a LaO and SrO layer. 
However, this increases the magnetic moment for Mn sandwiched between LaO by 0.04 iib only. A similar effect has 
been noted for calcium perovskitea^. 

VI. OXYGEN DISPLACEMENTS AND CHARGE ORDER IN LSMO 

We investigated the influence, on the Mn valency, of a tetragonal shift of the oxygen atom away from the SrO layer. 
An upward shift in the z-direction of 5% of the lattice constant was implemented. We calculated the total energies for 
a double unit cell, in the FM state, of La(i_3;)Sra;Mn03 (with a;=0.3 and 0.5) of the following configurations: all Mn^+, 
all Mn'^"'", and mixed Mn'*'*'/Mn^+ (Fig. Oil- For the mixed valence configuration, the Mn'^+ atom is taken to be inside 
the octahedron which is elongated by the tetragonal shift (see the cell on the right hand side of Fig. In the case of 
Lao.sSro.sMnOa, we found that increasing the tetragonal shift for the oxygen atom decreases the difference between 
the energy of the ground state (AF-A with Mn'*+) and the energy of the mixed valence Mn^+/Mn'*+ configuration. 

We have calculated a change in the AF Mn^+ ground state configuration as a function of tetragonal shift of O 
atom. For Lao.sSro.sMnOa the new mixed valence Mn^+/Mn^+ ground state occurs at 4% shift in the rigid band 
model and a 3% shift in the supercell model. In both models the Mn'^+ configured atom was placed in the elongated 
MnOg octahedron. Placing the Mn^+ in the elongated octahedron and the Mn'^+ in the squashed octahedron was 
the most unfavourable scenario. The crossover to a new ground state also depends on the Sr concentration and for 
Lao.ySro.aMnOa we find the new ground state already at 1% upward shift of O. 

The difference between rigid band model and supercell results for Lao.sSro.sMnOa shows the important role of the 
local environment, such as presence of Sr atoms, for the electronic properties of this material. Also, the magnetic 
moment on Mn'^^ increases slightly from 3.52 fiB to 3.63 /is, which is connected with O.le larger charge transfer from 
this atom, and decreases from 2.91 /is to 2.73 /is on Mn'*+. The oxygen atom, shifted in the supercell calculation, 
has a magnetic moment of 0.22 /xs, twice larger than in the undistorted structure. 
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In the bottom panel of Fig. |3]we present the polarisations of the electrons at the Fermi level. For shifts larger than 
3% of the lattice constant the system becomes fully polarised at the Fermi level and we obtain a half-metallic state. 
In Fig. ^we show the density of states for mixed valence Mn*+/Mn^+ system without tetragonal shift, mixed valence 
Mn''+/Mn'^+ system with tetragonal shift for oxygen atom (4% of lattice constant) and for AF-A system with Mn^+ 
configuration. When increasing the shift, the bottom of the conduction band moves up in energy to 0.07 eV above 
the Fermi level. This system then becomes half-metallic with a band gap in one spin-channel equal to 0.48 eV. 

In Fig. |31 we present a different tetragonal displacement involving two symmetrical shifts of 5% of the lattice 
parameter for two oxygen atoms from LaO and SrO layers in La2SrMn309. For the mixed valence ground state 
configuration Mn'^+/Mn^+/Mn'^+ we obtain a half-metallic state. In comparison with the above Mn*+/Mn'^+ double 
unit cell, the band gap has become nearly twice larger, at 0.8 eV. This is due to the higher concentration of Mn'^+ 
atoms in the latter supercell. The results of these two implementations of tetragonal distortions show that if random 
Jahn- Teller distortions occur, the system can become locally half-metallic. 

VII. PHASE SEPARATION AND CHARGE ORDER IN LSMO 

Supercells in the FM regime were also constructed to model both the influence of Mn''+ and Mn'^+ ordering on 
the total energy and the influence of the ordering of the LaO and SrO planes on the Mn charge order. Specifically, 
supercells of the form La„SrMn(„_|_i)03(„-|_i) with n = 2, ..,7 and La4Sr2Mn60i8, and La4Sr4Mn8024 were studied. 
Calculations for these systems were in agreement with the results of Fig. ^ Specifically, for n smaller than 4, the 
La„SrMn(„^i)03(„+i) systems had an Mn'*+ ground state and for n larger than 4, the mixed valence Mn^+/Mn'^+ 
system became the ground state. 

For the mixed valence configuration, we investigated the influence of different distributions of Mn*+ and Mn'^+ on 
the total energy. In particular, we studied different scenarios of distributions of Mn^+/Mn^+ atoms from all Mn^+ 
(denoted by (444444) in Table |llll) to aU Mn^+ (referred to as (333333) in Table HID, using both a six chemical units 
rigid band model of Lao.ssSro.iyMnOa and the LasSrMugOig supercell. In the supercell, the SrO layer was taken to 
be at the bottom of the cell. The energy differences are small and to emphasize this we write them down in meV. 
They are nearly equal to the magnitude of magnetic structure energy differences. The distribution 434343 of Mn''+ 
and Mn'^+ valencies was the most favourable energy state at 17% Sr concentration for both the supercell and the rigid 
band model calculations. A different distribution of the same amount of Mn''+ and Mn'^+ atoms denoted by 444333 
is about as unfavourable as localizing an extra Mn d electron, as realized in 334433 scenario, or delocalizing an Mn 
d electron as it happens in 443344 scenario. The energy differences between rigid band model calculations and the 
supercell show that SrO and LaO layers do matter. In the rigid band model we see smaller energy differences upon 
increasing the amount of Mn^+ cations (only 7 meV separate the 444444 state from the ground state) than upon 
decreasing the amount of Mn'^^ cations (215 meV separate the (333333) state from the ground state). In the supercell 
approach this trend is not as pronounced and whilst 213 meV separate the 333333 state from the ground state, the 
444444 configuration is further away from the ground state, by 27 meV, than in the rigid band model. 

The La4Sr2Mn60i8 system, i.e. increasing the Sr concentration to ~30%, acquires the Mn'*+ ground state config- 
uration and is metallic. For this system we have studied the influence of the distribution of SrO and LaO layers on 
the total energy. The state with the lowest energy was attained by phase separating the SrO layers from the LaO 
layers, namely the two SrO layers were nearest neighbors. Separating the two SrO layers by one and two LaO layers, 
respectively, increases respectively the energy of the system by 2 and 5 mRy. 

In Fig. IHl we show density of states (DOS) of one of the configurations of La4Sr2Mn60i8. The configuration is 
such that we have one Mn02 plane with Mn^+, sandwiched between two SrO planes, and all other Mn02 planes are 
occupied by Mn^+. This supercell can be considered as a model of phase separated LSMO which was studied by 
KoidaiSi It is also of relevance to the Lao.rCao.sMnOa system investigated by Bibesi^ As can be seen from Fig. 
this system LSMO is half-metallic. Energetically unfavourable, however, the energy difference between this state and 
the ground state is reduced to 20 mRy, in comparison with the energy difference between all Mn having the Mn'^^ 
configuration and the ground state (33 mRy). The gap is reduced from 0.68 eV (the scenario with the configuration of 
Mn'^+ for all manganese atoms) to 0.54 eV. Note that this pseudogap is constant for all Mn02 layers, independently 
of their valence. None of the supercells had this Mn'^^ valence in their ground state configuration. However one could 
speculate that maybe at the surface Mn'^^ could be the stable configuration since the lower coordination would favour 
the more localised state. More complex supercells, where more Sr layers could phase separate, could plausibly lead 
to a half-metallic ground state. 

The magnetic moment increases from 2.82 /is for Mn sandwiched between two SrO planes (denoted by MnsrO-Sro) 
to 2.99iiB for MuLaO-LaO (Mn sandwiched between two LaO planes) in the Mn'*+ configured ground state and from 
3.45/iB for MnsrO-SrO to 3.67 fiB for MuLaO-LaO in the half-metallic Mn'^+ configuration. We note that the magnetic 
moments on the oxygen atoms in the SrO layers are dramatically reduced in the vicinity of Mn''+ , namely, a value 
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TABLE III: Total energy differences (in meV), with respect to the state with the lowest energy, for different distributions of 
manganese ions, in six units supercell in the rigid band (corresponding to Lao.gsSro.irMnOs) and, in the superceU LasSrMnsOig. 
For LasSrMngOig the system is ordered as SrO-LaO-LaO-LaO-LaO-LaO and the Mn*"*" or Mn^^ ions are sandwiched in between. 

Scenario Supercell Rigid Band Model 
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of 0.01/iB is obtained, and they have values of 0T5 and 0.2/iB in the vicinity of Mn'^+/Mn'*+ and Mn^+ respectively. 
The magnetic moments of the O atoms residing in the LaO plane remain small at 0. OS/is independently of the Mn 
valency. In contrast, the O magnetic moments in the Mn02 planes are small with values around 0.05/iB and do not 
change much depending on Mn valency except for Mn^+ at the SrO/LaO interface were an O magnetic moment as 
large as O.l/xs has been calculated. 

Through the study of the La4Sr4Mn8024 supercell we model an even further increase of the Sr concentration 
to 50%. The geometry of the supercell is four La layers separated from the four Sr layers. For this supercell 
the ground state configuration is all Mn^"^ and this has 31 mRy lower energy than the half- metallic system with 
mixed valence Mn^+/Mn'^+ configuration (33333444) and 61 mRy lower than the configuration with all Mn'^+. The 
magnetic moments, which are smaller than for 30% Sr, are increasing, as is also the case for 30% Sr, from 2.71 fiB for 
MnsrO-SrO to 3.05 hb for MuLaO-LaO in the Mn''+ configured ground state, and from 3.56/iB for MnsrO-SrO to 
3.61/zb for MuLaO-LaO in the half-metallic Mn'^+ configuration. For the latter, more localised configuration, we note 
a large sensitivity of the magnetic moments of the oxygen atoms from the LaO layers to the presence of the SrO /LaO 
interface. At the interface the O magnetic moment of the LaO layer practically disappears. Away from the interface, 
this moment is 0.1 /is and remains constant in these LaO layers. In the SrO layers the O magnetic moment is slightly 
larger, at O.lS/is, and remains constant in all SrO layers including the interface one. 

VIII. CONCLUSIONS 

We have found that the ground state configuration of LSMO consists of manganese in the 4-1- valence configuration. 
For Sr concentrations less than 20%, a mixed valence Mn'*+/Mn^+ ground state was obtained, whilst for concentrations 
in excess of 20% a Mn'*+ ground state was obtained. Also, a close competition between FM and AF-A Mn"'+ states 
was seen with a crossover from FM to AF-A Mn^+ around 35% Sr concentration. The Mn'^^ ground state has 
marginal spin polarization at the Fermi level. We found, however, that the Mn'^^ valency was not that energetically 
unfavourable and that small increases in the lattice parameter would lead to the occurence of Mn'^^ in a mixed valence 
Mn''+/Mn'^+ state. Likewise, elongations of the MnOe octahedra lead to the formation of Mn'^^ valence atoms. With 
the appearance of Mn'^"'" ions a half-metallic state is obtained. Furthermore, in the mixed valence state we found that 
the total energy is minimized by an ordered array of Mn^+ and Mn^+ Mn02 planes. The relation of this finding to 
the observed charge-ordered stripes^ will be the topic of a further study. 

As the Sr doping increases, the Mn''+ state becomes more and more favoured with respect to the Mn^+. The Sr 
doping results in one more valence band electron which is obtained from delocalizing a Mn'^+ Cg state. In other words, 
Sr hole doping favours band formation instead of localization. 

Finally, for Lao.ySro.sMnOs, of huge relevance to spintronics applications, we suggest the following scenario to 
explain the half-metallic properties. Whilst in Lao.rSro.sMnOa, from x-ray results, every Mn-0 bondlength is the 
same at 1.949 A,^^ the bondlengths of Mn-0 in distorted LMO change from 1.907 A to 2.179 A, depending on 
direction. It seems however, based on pulsed-neutron diffraction measurements that the Jahn- Teller distortion is 
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still present in LSMO and could therefore, according to our analysis, lead to the formation of Mn^+ ions, which would 
be sufficient to give rise to a half-metallic state. 
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FIG. 1: Left: Energy phase diagram LMO — LSMO in the cubic structure with the lattice constant of 7.32 a.u.. The energy 
for all the systems is with respect to the energy of FM system with Mn'*"''. Top: AF-A — crosses with solid and dashed lines 
for Mn^"*" and Mn*"'" respectively; AF-G — squares with solid and dashed lines for Mn^"*" and Mn''"'" , respectively. Bottom, FM: 
Mn^"*" - squares with solid line; Mn*"'"/Mn^"'" supercell — solid circles with solid line; Mn^+ZMn^^/Mn^"*" supercell - triangles 
with solid line; Mn*"'"/Mn^"'"/Mn''"'" supercell — triangles with dashed line; Mn^"*" for the lattice constant of 7.434 — squares 
with dashed line; Mn^+ for appropriate supercell — open circles with solid line. The crosses on their own show AF-A ground 
state for LSMO. Right: Schematic views of the AF-A and AF-G structures. 
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FIG. 2: Density of states for FM La.rSr.sMnOa from rigid band model for minority and majority spin channels with respect 
to the Fermi energy. Displayed are: LSD calculation (top), SIC— LSD calculation for the ground state of three localised t2g 
electrons (centre) and SIC— LSD calculation for four localised electrons (3t29+eg(3^2_^2)) (bottom). 
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FIG. 3: Left: In the top panel the energy difference for the FM system in the Mn*"'"/Mn^"'" mixed valence state: squares with 
solid line - rigid band model (for x=0.5); open circles - supcrcell model (LaSrMn206); squares with dotted line - rigid band 
model (for x=0.3), and AF-A with Mn''+ - crosses with dashed line. All energies arc with respect to the energy of the FM 
system with Mn'*+. In the bottom panel: polarisation of electrons at the Fermi level: squares - from rigid band model (for 
x=0.5) and open circles for supercell model (LaSrMn206). Right: Schematic view of the unit cell. 
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FIG. 4: DOS for LaSrMn206 with mixed valence Mn'''*"/Mn''^ configurations for minority and majority spin channel: is 
metallic without tetragonal shift for oxygen atom (top), is half- metallic with a tetragonal shift applied to the oxygen atom of 
0.04 of the lattice constant (centre). In the bottom panel is the DOS for the antiferromagnetic ground state of Mn*"*" without 
tetragonal shift. 
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FIG. 5: Left: DOS, in minority and majority spin channels, for La2SrMn309 with symmetrically shifted (by 0.05 of lattice 
constant) oxygen atoms for the mixed Mn^"'"/Mn'*"'"/Mn^"'" configuration. Right: Schematic view of the unit cell. 
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FIG. 6: Local DOS for a supcrccll La4Sr2Mn60i8 (right). The configuration is (a) one Mn02 layer of Mn*^ sandwiched 
between the SrO layers (marked as black balls) and Mn^+ in all other Mn02 planes. Here (b) refers to SrO plane, (c) to Mn02 
plane, and (d) to LaO plane. The left-hand-side picture shows the structure. 



